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ABSTRACT 

We report a new technique to select 1.6<z<4.6 dusty Lyman-alpha emitters (LAEs), over a third 
of which are 'blobs' (LABs) with emission extended on scales of 30-100 kpc. Combining data from 
the NASA Wide-field Infrared Survey Explorer (WISE) mission with optical spectroscopy from the 
W. M. Keck telescope, we present a color criteria that yields a 78% success rate in identifying rare, 
dusty LAEs of which at least 37% are LABs. The objects have a surface density of only ~ 0.1 deg -2 , 
making them rare enough that they have been largely missed in narrow surveys. We measured 
spectroscopic rcdshifts for 92 of these WISE-selected, typically radio-quiet galaxies and find that the 
LAEs (LABs) have a median rcdshift of 2.3 (2.5). The WISE photometry coupled with data from 
Hcrschcl 11 reveals that these galaxies have extreme far-infrared luminosities (Ltr . > 10 13-14 Lo) and 
warm colors, typically larger than submillimeter-selected galaxies (SMGs) and dust-obscured galaxies 
(DOGs). These traits are commonly associated with the dust being energized by intense AGN activity. 
We hypothesize that the combination of spatially extended Ly-a, large amounts of warm IR-luminous 
dust, and rarity (implying a short-lived phase) can be explained if the galaxies are undergoing strong 
'feedback' transforming them from an extreme dusty starburst to a QSO. 

Subject headings: galaxies: high-redshift — galaxies: starburst — infrared: galaxies — galaxies: ISM — 
galaxies: formation 



1. INTRODUCTION 

High-redshift Ly-a emission is widely used to study 
star formation in galaxies (e.g. Cowie & Hu 1998; 
Stcidcl et al. 2000). Systems that exhibit this line 
are typically referred to as Ly-a emitters (LAEs); a 
small subset show extended emission on scales > 30 kpc 
and are considered Ly-a 'blobs' (LABs). Among the 
largest coherent galactic structures known in the Uni- 
verse, LABs are extremely energetic (with Ly-a lu- 
minosities of ~T0 42_44 ergss _1 ) and have been shown 
to trace over-densities of galaxies at high-redshift (e.g. 
Keel et al. 1999; Francis et al. 2001; Steidel et al. 2000; 
Matsuda et al. 2004; Steidel et al. 2011). Optically- 
selected LABs are 100-1000 times less abundant than 
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LAEs with only a few dozen known (e.g. Stcidcl et al. 
2000; Matsuda et al. 2004). 

Resonant scattering by atomic hydrogen results 
in Ly-a photons being easily extinguished by dust. 
Hence, blind narrow-band searches for z ~ 2 — 6 LAEs 
and LABs have been based on the paradigm that 
high-redshift LAEs are generally blue galaxies, with 
little to no dust (Ay=0.0-2.5 mag), and moderate star 
formation rates (10-100 M Q yr _1 ; e.g. Cowie & Hu 1998; 
Rhoads et al. 2000; Steidel et al. 2000; Matsuda et al. 
2004; Gawiser et al. 2007; Gronwall et al. 2007; 
Nilsson et al. 2007, 2011; Finkelstein et al. 2009). 

This suggests that dusty infrared galaxies and 
LAE/LAB populations, despite both peaking in activ- 
ity around z ~ 2, should largely not overlap. In- 
deed Spitzcr observations of ~ 40 optically selected 
z = 2 — 3.5 LABs found that < 15% are luminous in 
the infrared (Webb et al. 2009; Nilsson & M0ller 2009; 
Lai et al. 2007; Geach et al. 2005). Furthermore, none 
of the optically discovered LABs (e.g. Steidel et al. 2000; 
Matsuda et al. 2004, 2011) are detected at 12 or 22/um 
by WISE. However, 40-50% of the infrared luminous 
sub-millimctcr galaxies (SMGs) show Ly-a emission, al- 
though with typically moderate line fluxes and rarely 
with spatially-extended Ly-a emission (Chapman et al. 
2005). This suggests that dust and Ly-a emission are 
not in fact mutually exclusive, and points to a clumpy 
dust distribution that allows sight lines for Ly-a pho- 
tons to escape (Ncufcld 1991). Since SMGs and related 
high-redshift galaxies are thought to be just one stage 
in the evolution of massive elliptical galaxies, the study 
of dusty galaxies with spatially extended Ly-a can thus 
provide a unique insight to this process. 
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This paper, the first in a series, presents a new, highly 
efficient mid-IR color technique to select a unique pop- 
ulation of z — 2 — 4 dusty LAEs and LABs. We use 
the redshift distribution, general spectroscopic proper- 
ties, and preliminary Herschel observations to deter- 
mine how these objects relate to other dusty galaxies 
at these redshifts. We assume a ACDM cosmology with 
Q M = 0.27, fl A = 0.73, and H =71 kms -1 Mpc -1 and 
use Vega magnitudes. 

2. SAMPLE DEFINITION AND OBSERVATIONS 

One of the primary science goals of the NASA 
Wide-field Infrared Survey Explorer (WISE) mission 
(Wright et al. 2010) is to find the most luminous galax- 
ies in the Universe. Initial searches in WISE color space 
revealed a diversity of galaxy types over a wide range 
of redshifts (Griffith et al. 201 1; Wu et al. submitted). 
Deep spectroscopic follow-up of the galaxies with the red- 
dest mid-IR colors often showed intense Ly-a emission 
(also see Eiscnhardt ct al. (submitted)). A striking fea- 
ture of these WISE-sclccted LAEs, is that at least one- 
third are blobs (LABs), having spatially extended Ly-a 
reaching 30-100 kpc (see Figure I). 

2.1. Selection Method 

Realizing that WISE-detected LAEs and LABs lie in a 
very specific region of a WISE color-magnitude diagram, 
we refined the selection criteria used in the initial extra- 
galactic follow-up to focus on high- redshift LAEs/LABs. 

The resulting selection criterion, shown in Figure 2, 
are W2-W3 (4.5-12^m)>4.8 and SNR >5 in W3 and 
W4. We also excluded sources within 30 degrees of the 
Galactic center and 10 degrees of the Galactic plane to 
avoid AGB stars and saturation artifacts. In order to 
remove low-redshift (z 1 0.5) star-forming galaxies with 
similar WISE colors, we required non-detections in the 
SDSS and DSS imaging (r'>22). Finally, each object's 
coadded W3 and W4 images were visually inspected to 
ensure the source is not spurious. These selection criteria 
result in —4400 galaxies over the whole sky. 

Another WISE color selection, "W12drop", outlined 
in Eisenhardt et al. (submitted) has been successful in 
identifying ULIRGs over a wider redshift range (z ~ 
0.05 — 4.6), and while there is substantial overlap in the 
2> 2 samples, overall the refined WISE-LAE/LAB crite- 
rion is at least 2 (2.5) times more effective at selecting 
high-z LAEs (LABs). The primary differences arc that 
the WISE-LAB criterion discussed here impose an opti- 
cal magnitude limit to remove contamination from low 
redshift galaxies, place no flux cut on Wl, and probes a 
factor of 1.2-2 deeper in W3 and W4 and 0.5 mags bluer 
in W2-W3 than the W12drop selection. 

2.2. Optical Spectroscopy 

To determine redshifts and discern the nature of the 
galaxies, we used Keck I LRIS (Oke et al. 1995) over 
the course of six runs between July 2010 and Jan- 
uary 2012. The sample contains 101 objects that fulfill 
the WISE LAE/LAB selection, of which -45% had al- 
ready been followed-up using the W12drop criteria, and 
therefore did not need to be re-observed. The major- 
ity of data were taken using the 600 £ mm -1 grism in 
the blue arm (Abiazc =4000A; spectral resolving power 



R = A/AA -750), the 400^ mm" 1 grating on the red 
arm (A b iazo =7800 A; R -700), the 5600A dichroic and a 
l'/5 wide longslit. The targets were observed with a me- 
dian seeing of Of! 7 and total integration times between 
20-40 minutes. The LRIS data were reduced and flux 
calibrated using standard procedures. 

The optical spectroscopy reveal a wide range of char- 
acteristics, from pure starburst-dominated galaxies with 
narrow emission lines and several ultraviolet interstel- 
lar absorption lines comparable to Lyman break galaxies 
(Stcidel et al. 2000) to those with strong AGN compo- 
nents, inferred from the presence of high ionization lines 
such as NV, SilV/OIV, CIV and broad emission line pro- 
files of several thousand kms -1 . 

Robust redshifts based on two or more spectral fea- 
tures were determined for 92 out of the 101 sources we 
attempted (91% success rate), with a redshift range of 
z = 1.13 — 4.59 (see Figure 3), while the remaining 9 
showed no continuum or spectral lines. The excellent 
blue sensitivity of the LRIS detector allows Ly-a emis- 
sion to be detected down to 3160A, corresponding to a 
redshift of z —1.6. Of the 92 WISE selected galaxies with 
a robust redshift, 89 are at a z > 1.6, and of those 79 
(89%) showed Ly-a in emission. An additionally strik- 
ing finding was that 37% (29/79) of the WISE LAEs 
showed Ly-a emission in the two-dimensional spectra 
extended on spatial scales of > 30 kpc (i.e., LABs), with 
18% (14) having emission 40-100 kpc in extent. We also 
note that 14 have extended emission of 25-30 kpc, which 
is still more than 5 times the size of the galaxies, and 
more than twice the seeing limitations, constituting small 
LABs or Ly-a clouds. 

The definition of a LAB does vary between stud- 
ies. Matsuda et al. (2004) required an isophotal 
area of at least 16arcsec 2 (or —30 kpc in length 
at z = 3.1) and a surface brightness of —27.6 
mag arcsec -2 (— 2 x 10 -18 ergss -1 cm" 2 arcsec -2 ), 
while Matsuda et al. (2011) probed down to 
1.4xl0 -18 crgss -1 cm -2 arcsec -2 . Making basic as- 
sumptions (slit width and size of extended emission) we 
measure a surface brightness of the > 30 kpc-scale blobs 
of 1-100 xlO -18 ergs s -1 cm -2 arcsec -2 which is slightly 
brighter than the limits used in the optical surveys. It 
is therefore plausible that the true sizes of the WISE 
LABs are in fact larger under these standard definitions. 

2.3. Far-Infrared Photometry 

The red WISE colors strongly argue that significant 
amounts of warm dust are present in these galaxies. It 
is common to compare z — 2 dusty galaxies by charac- 
terizing their dust luminosity and temperature. Since 
their spectral energy distributions (SEDs) typically peak 
in the far-infrared, observations between 0.1 — 1 mm are 
required. To that end, we have undertaken a campaign 
with Herschel to study the WISE-selected LABs in five 
bands using the PACS (Poglitsch et al. 2010) and SPIRE 
(Griffin et al. 2010) imagers. The PACS data were ob- 
tained using the mini scan-map mode, and reached a 
5cr depth of 30mJy at 160^m. SPIRE observations 
employed the small jiggle-map mode, and reached the 
250/im confusion limit of — 8mJy. The data were re- 
duced and photometry extracted using recommended 
procedures with the HIPE software (Ott 2010). 
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Fig. 1. — Two-dimensional LRIS spectra for four of the WISE 
LABs in the sample. The spatial extent of the Ly-a emission is 
noted in the corner of each image. These examples are representa- 
tive of the asymmetric, spatially extended line profile, and clumpy 
nature of the Ly— a emission found in these galaxies. 

Of the 8 objects observed so far 2 , all were detected at 
160/im and of those, 6 are also detected at 250/im. We 
fit these data using a simple SED based on a modified 
blackbody and in Figure 4 compare far-IR luminosity and 
dust temperature of the WISE LABs against SMGs and 
DOGs from the catalogs of Magnelli ct al. (2012) and 
Melbourne et al. (2012). It is clear that the WISE LABs 
typically contain warmer dust, suggesting that they are 
being heated, at least in part, by an AGN, since star- 
burst galaxies such as SMGs are typically colder (also 
see Wu et al. (submitted)). The WISE LABs are sys- 
tematically more luminous than SMGs, with bolometric 
luminosities in excess of Lpm > 10 13 L Q , making them 
Hyper-Luminous Infrared Galaxies (HLIRGs). 

It might be possible to associate these high luminosities 
with gravitational lensing; however, none of the optical 
spectra demonstrate contamination from a foreground 
source, and the spatially extended Ly-a emission for all 
the sources is at the same redshift as the other spectral 
features. Furthermore, adaptive-optics near-IR imaging 
of eight LABs with Keek's NIRC2 camera (Bridge et al. 
in prep.; Eiscnhardt et al. submitted) reveal no poten- 
tial lens (nearby galaxies or clusters) or shearing of the 
source. 

3. DISCUSSION 

One of the main features that sets the WISE LAEs 
apart from other high-redshift galaxies is the unprece- 
dented high fraction that exhibit spatially extended Ly-a 

2 A detailed treatment of the Hcrschcl properties of WISE LABs 
is deferred to a separate paper when more data is collected. 
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Fig. 2.— WISE W3(12/xm ) vs. WISE W2-W3 (4.5-12/mi) color 
for all sources in the sample with spectroscopic redshifts between 
1.6 < z < 4.5. WISE LAEs (black circles), 30-40 kpc LABs (green 
squares), and 40-100kpc LABs (orange squares) are plotted. For 
comparison, z = 2 — 3 WISE detected SMGs and DOGs arc high- 
lighted. The WISE LAEs (LABs) exhibit unique, redder mid-IR 
colors (implying hotter dust temperatures) than other dusty z ~ 2 
populations. W2-W3 color selection criteria has a 78% (>29%) suc- 
cess rate in identifying z > 2 dusty LAEs (LABs) (dashed line). 
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Fig. 3. — Spectroscopic redshift distribution for the full sam- 
ple (LAEs+LABs; filled grey), and WISE LABs (filled orange). 
For comparison we plot the redshift distribution of other well- 
studied dusty populations: SMGs (black dashed; Chapman et al. 
(2005)), Herschel SMGs (blue long dashed; Magnelli et al. (2012)), 
and DOGs (red dotted; Bussmann et al. (2012)). 



emission, despite having large amounts of dust. We cross 
matched the WISE catalog with all previously known 
optically selected LABs in the literature and found that 
none are detected at 12 or 22/im above a SNR of 5, in- 
cluding the ones known to host ULIRGs (Geach et al. 
2005; Prescott et al. 2009). Despite being a factor of 10- 
1000 times more luminous in the mid-far-IR, the WISE 
LABS have Ly-a luminosities (10 42 ~ 44 ergss -1 ), rest- 
frame equivalent widths (6-300A), and blobs sizes simi- 
lar to optically discovered LABs. With a surface density 
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Fig. 4. — Estimated dust properties based on Herschcl observa- 
tions. Objects between 2.0 < z < 3.5 extracted from catalogs of 
DOGs and unlcnscd SMGs are plotted against the eight observed 
WISE LABs over the same redshift range. In all cases, fits are to 
a single temperature modified blackbody. WISE objects without a 
detection at 250/^m are plotted as open circles, and represent lower 
limits to both the temperature and luminosity. Error bars in all 
cases are based on fit results only and do not account for other 
systematic effects. 

of 0.03 per square degree however, the WISE LABs are 
thousands of times more rare. 

3.1. WISE LABs and Other Dusty High- Redshift 
Populations 

Though the IR luminosity and redshift distribution of 
the WISE LAEs/LABs are similar to other populations 
of IR- luminous galaxies (Figure 3), they must differ in 
some way given the large fraction with extended Ly-a 
emission. We now compare the WISE LAEs/LABs to 
other high-redshift dusty populations and also note how 
the various selections overlap. A more in-depth com- 
parison will be presented in the upcoming catalog paper 
(Bridge ct al. in prep.). 

Blind sub-millimeter surveys have been very successful 
at finding high- redshift, dust-obscured galaxies with far- 
IR luminosities in excess of 10 12 L Q (i.e. SMGs). How- 
ever, these long wavelength observations (^500/im) se- 
lect systems with cooler dust and less extreme mid-IR 
colors than seen in our WISE sample (Figure 2 & 4). 
Generally, SMGs are powered by star formation rather 
than Compton-thin AGN (Alexander et al. 2005). Based 
on the limited Herschel observations conducted so far, all 
the WISE LABs would also be classified as SMGs (Fig- 
ure 4). However only 1-2% of known SMGs fulfil the 
WISE LAE/LAB criteria, and there are only a few rare 
SMGs with extended Ly-a emission. 

Dey et al. (2008) have successfully used a mid-IR to 
optical flux ratio {FzAtxm/FR >1000 and i^/xm > 
0.3mJy) to select a population of z ~ 2 dust-obscured 
galaxies (DOGS; Houck et al. 2005; Desai et al. 2009; 
Bussmann et al. 2012). By probing blueward of the 
SED peak, this technique recovers a mix of AGN and 
starburst dominated galaxies. Although all the WISE 
LAEs/LABs fulfil this criteria (using r'-W4), only < 8% 
of the DOGs in Bussmann ct al. (2012) fulfill the WISE 
selection presented here. WISE LAEs/LABs are also sig- 
nificantly brighter in the mid-far-IR (>5 times). Finally, 



the DOGs criterion does not specifically select sources 
with a steeply rising mid-IR SED (Fig. 2), that we have 
shown are generally the objects that exhibit extreme 
Ly-a properties. 

Our sources are typically undetected or only faintly 
detected in the FIRST (<lmJy) or NVSS radio surveys, 
differentiating them from the well-studied high-redshift 
radio galaxies (Reuland et al. 2003) and QSOs which can 
also exhibit large Ly-a halos (e.g. McCarthy et al. 1987; 
Smith et al. 2008). 

In short, though surveys of DOGs/SMGs will also find 
WISE LAEs/LABs, with a projected source density of 
roughly 1 per 30 deg 2 , only large-area/all-sky surveys 
are able to find a statistically significant sample. 

3.2. Spectral Properties of WISE LABs 

The cause of the LAB phenomena in optically se- 
lected samples is still debated (e.g. Prescott et al. 2009; 
Colbert et al. 2011), but given that all WISE LABs 
are associated with a luminous IR galaxy while less 
than 15% of optical LABs arc (e.g. Webb et al. 2009; 
Nilsson & M0llcr 2009), the dominant powering mecha- 
nisms behind these two species of LABs are likely differ- 
ent. For clues we investigate the optical spectra. 

The high-ionization and broad emission lines seen in 
Z 70% of optical spectra suggest that the majority, if 
not all, of the WISE LAEs/LABs host a dust-obscured 
AGN, which is also consistent with the Herschel and 
WISE colors. The Ly-a emission is spatially asymmet- 
ric and clumpy, and is often offset from the systemic 
redshift of the galaxy, implying large-scale flows. Fur- 
thermore, many galaxies in the sample demonstrate a 
Ly-a velocity structure rcdshiftcd up to several thousand 
kms -1 as a function of distance from the central compo- 
nent (Figure Id). The line profiles are diverse, ranging 
from traditional P-Cygni profiles characteristic of Wolf- 
Rayet and O-star winds to Ly-a emission peaking on the 
blue side of the Ly-a profile, suggesting possible inflows 
(Dijkstra & Loeb 2009; Barnes et al. 2011). 

Multiple slit orientations were used to probe the spatial 
line morphology for four of the LABs. A slit position 
angle rotated by 45 and 90 degrees of the original angle 
often showed little or no extended emission, implying an 
asymmetric and filamentary line morphology. This result 
suggests that the quoted efficiency of the color selection 
in identifying LABs is a lower limit. Indeed, using a 
simple model of the extended emission size and geometry, 
and assuming that all LAEs in our sample are LABs, 
the 37% detection rate is consistent with a random slit 
orientation. 

3.3. WISE LABs: Feedback Caught in the Act? 

The correlation between the masses of supcrmas- 
sive black holes (SMBH) at the center of nearby el- 
lipticals and their bulge stellar velocity dispersion re- 
mains striking in the study of galaxy evolution (e.g. 
Ferrarese & Merritt 2000). The general paradigm that 
ties these passive galaxies with their high-redshift pro- 
genitors is a process of merger-induced star formation 
which also fuels the SMBH (e.g. Sanders & Mirabel 1996; 
Hopkins ct al. 2006). The system appears at one point 
as a heavily obscured star forming galaxy, and later, 
as the SMBH accretes at a higher rate, a galaxy with 
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an active nucleus. AGN- and starburst-induced winds 
eventually become strong enough to expel the obscuring 
gas and dust, briefly revealing an optical quasar. This 
short-lived 'feedback' process is thought to quench both 
star formation and AGN activity, leading ultimately to 
a passive, red galaxy spheroid (e.g. Hopkins et al. 2006; 
Farrah et al. 2012). Observational evidence for this pro- 
cess has been challenging to obtain, but models provide 
a general picture of the properties of these systems. 

It is well established that the peak of this feedback 
activity occurs near z ~ 2 (e.g. Di Matteo et al. 2005), 
and that a multitude of galaxy types are associated with 
the process. SMGs are thought to be the IR- luminous 
starbursting prc-cursors to systems that will evolve into 
a massive elliptical, and have been well characterized in a 
number of studies (e.g. Chapman et al. 2005; Borys et al. 
2005; Magnelli et al. 2012). 

Using a 3D radiative equilibrium code, 
Chakrabarti et al. (2007) showed that AGN feed- 
back is particularly effective at dispersing gas and dust. 
They predict that for a brief time (< 40Myr) while the 
AGN injects energy into the galaxy and surrounding 
IGM at its maximal rate, the dust is hotter than during 
the SMG starburst phase. Strikingly, all the WISE 
LABs observed with Herschel thus far show these warm 
far-IR colors consistent with being in a short-lived AGN 
feedback phase. The outflow of material would not 
necessarily be symmetrical, and Ly-a could appear 
filamentary. 

The WISE LABs which demonstrate warm, IR- 
luminous dust and directional Ly-a emission, share 
these general properties. The surface density of LABs is 
generally consistent with predicted numbers of HLIRGs 
from models (e.g. Bethermin et al. 2011) coupled with 
timcscalc arguments for the feedback process and length 
of the SMG phase. However we caution that this is not 
a particularly strong argument given the uncertainty in 
the models. Nevertheless, the accumulated evidence sug- 
gesting that the WISE LABs are undergoing extremely 
powerful feedback is intriguing, and we will expand on it 
in future papers. 

4. CONCLUSIONS AND FUTURE WORK 

This letter presents a new WISE color criterion with a 
78% success rate in selecting z= 1.6 — 4.6 dusty, mid-IR 
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bright, Ly-a emitters, of which at least 37%, arc blobs. 
This new population of galaxies, largely missed in nar- 
rower surveys, have a rcdshift distribution that peaks at 
z ^2.3, total IR luminosities on average brighter than 
SMGs and other high-redshift dusty galaxies, and a den- 
sity of ~ 0.1 deg -2 . This is the first systematic search 
technique to highlight dusty LAEs/LABs, and unlike op- 
tical narrow-band searches covers a large redshift range, 
the whole sky, and without contamination by [Oil] in- 
terlopers, making it well suited to providing targets for 
future large spectroscopic surveys. 

We speculate that these galaxies are in a short-lived 
transition between a dusty starburst and an optical QSO 
driven by the central AGN. If true, these systems of- 
fer a unique opportunity to investigate AGN feedback 
and how it can effect not only the galaxy but also 
the surrounding intcrgalactic medium. A full census of 
the WISE LAEs/LABs, optical spectroscopy, mid-far-IR 
properties, and near-IR morphologies will be presented 
in a series of forthcoming papers. 
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